The origin of increase of damping in transition metals with rare earth impurities 
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Abstract 

The damping due to rare earth impurities in transition metals is discussed in the low concentration limit. It is shown that 
the increase in damping is mainly due to the coupling of the orbital moments of the rare earth impurities and the conduction 
p-electrons. It is shown that an itinerant picture for the host transition ions is needed to reproduce the observed dependence 
of the damping on the total angular moment of the rare earths. 

PACS numbers: 72.25.Rb, 76.30.Kg, 76.60.Es 
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Magnetization dynamics has become one of the most 
important issues of modern magnetism. This develop- 
ment is driven by the technological demand to tailor 
magnetic responses on ever smaller length and shorter 
time scales. The importance of this issue manifests it- 
self in a completely new area of research, spintronics, 
and a huge literature that cannot be cited here. Selected 
highlights include precessional switching by tailored field 
pulses 0, @ , spin-torque 0, 0] , and laser-induced magne- 
tization dynamics [a, la] . 

In general, magnetization dynamics is described via 
the Landau-Lifshitz-Gilbert equation (LLG) Q includ- 
ing additional terms to incorporate spin-torque effects 
U or those due to pulsed optical excitations . All 
these descriptions account for energy dissipation via a 
phenomenological damping parameter a which governs 
the time needed for a non-equilibrium magnetic state to 
return to equilibrium. Recently it has even been sug- 
gested that a determines the magnetic response to ultra- 
fast thermal agitations [lfj. 

Technological applications call for the ability to tai- 
lor a [lH . The most systematic experimental investiga- 
tion on this topic was published by Bailey et al. |l2l | who 
studied the effect of rare-earth doping on the damping 
in permalloy. Most rare earth ions induced a large in- 
crease of a, but neither Eu nor Gd altered the damping 
of permalloy (cf. Fig.0) ■ Since Gd +++ and Eu ++ have 
no orbital momentum, this points immediately to the im- 
portance of the angular momentum in the damping pro- 
cess. Bailey et al. determined damping by reproducing 
their data via the LLG equation using a as a fit parame- 
ter. This widely used procedure points to a fundamental 
problem of this phenomenological approach. Though the 
LLG equation describes data well, a more microscopic 
approach is needed to understand the origin of damping. 

It was Elliott fl3j who first studied damping in semi- 
conductors due to spin-orbit coupling. Later Kambersky 
[l4| argued that the Elliot- Yafet mechanism should be 
also operable in magnetic conductors. Korenman and 
Prange fl5| developed a more microscopic treatment and 
found that spin-orbit coupling should be important at 
low temperature in transition metals. Recent measure- 
ments of damping in magnetic multilayers at room tem- 



perature suggest that the s-d interaction might also 
be at the origin of damping 0, ^gj. Howevever, all of 
the present models fail to reproduce the data of Ref. 

In this Letter, we explain the increase of damping in 
rare earth doped transition metals via a novel orbit-orbit 
coupling between the conduction electrons and the im- 
purities. The well known s-f interaction |19| gives rise 
to a (gj — l) 2 dependence of the dampingthat is in con- 
tradiction to experimental observations |12j . In contrast, 
the orbit-orbit coupling considered here reproduces the 
measured (gj — 2) 4 dependence of the damping. Both de- 
pendencies on the Lande g-factor gj follow directly from 
the fact that the rare earth ions are in their ground state. 
Hence, their angular momentum Lf, spin Sf, and total 
angular momentum Jf are related by the Wigner Eckard 
theorem: Lf = (2 — <?,/)Jf and Sf = (gj — l)Jf. Deriv- 
ing the magnetic moments of the transition ions from the 
electronic degrees of freedom is essential to capture the 
correct behavior of damping as a function of J/ . For the 
uniform mode, the damping due to orbit-orbit coupling 
is of Gilbert form in the low frequency limit. 

Taking the wave functions of the d-, /-, and conduction 
electrons orthogonal, the Hamiltonian for the rare-earth 
doped transition metal in an external field H is 



TL = TL e + TLf + Ttd- 



(1) 



This approximation should be valid for the heavy rare 
earths but probably fails for elements like Cerium 
where valence fluctuations are important. The conduc- 
tion electron Hamiltonian Ti e is the usual one, TL e = 
J2k.o- e k,o-ak,Ja k ,o-, where a\ t<J and a k , G are the creation 
and annihilation operators of a conduction electron with 
momentum k and spin a. e k ,a is the energy of the con- 
duction electrons including a Zeeman term. 

Jif is the Kondo Hamiltonian |2(| of the localized rare 
earth moment 



S e /f and L, 



= rS e • S/ + AL e • L 



/ - W 



H. 



(2) 



are the spin and angular momentum 



of conduction and / electrons, respectively. ^ e /f are 
taken with respect to the position of the impurity. The 
spin-spin term is the well known s-f coupling used by 
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de Gennes to reproduce the Curie-temperatures in rare 
earths with T being of the order 0.1 eV |19|, The last 
term is again a Zeeman term. The middle term is 
the essential orbit-orbit interaction needed in our dis- 
cussion. To get a non-zero orbit-orbit term due to 
a single impurity at the center, it is essential to in- 
clude higher terms of the partial wave expansion for the 
wave functions of the conduction electrons: ifiki?) — 
% EZo Y2=-i i l mMkr)Y lm (0 k , <f>). The 

first non-trivial contribution for I = 1 is j20| 

Hll = i(2~ gj)J2\ (k, k')kxk'-J f a\a k , , (3) 

k,k' 

where the orbit-orbit coupling A will be assumed to be 
a function of the relative angles of the k vectors and is 
almost everywhere zero except for k close to the Fermi 
level hp. The magnitude of A is not known but is ex- 
pected to be of the same order as the spin-spin coupling 
constant T [2ll . The crystalline electric field effect in 
transition metals is less than 0.1 meV which is small and 
hence the spin-orbit term S e ■ is neglected. At room 
temperature all the rare earth ions studied in Ref . [12J are 
in their ground state making the term Sf ■ ineffective 
as damping mechanism. This follows immediately from 
the Wigner-Eckart theorem. 

The Hamiltonian for the host transition ions is based 
on the Anderson Hamiltonian with explicit spin ro- 
tational invariance in the absence of a Zeeman term 
[THHIH. It is 

Hd = e d d^d a + 2J V k d (ah,ad a + d\a k ^) 

k 

+ jP 2 - jSd ■ Sd - Hd • H, (4) 

where Sd is the spin operator of the local d electrons while 
their orbital angular momentum is assumed quenched. 
p is the charge density operator of the d electrons. In 
transition ions such as Ni, Vkd ~ 1-0 — 10.0 eV is com- 
parable to the Coulomb potential U. The hybridization 
term between the conduction- and d-electrons is essential 
to establish a spin-independent orbit-orbit coupling be- 
tween the d- and the f-ions. The degree of localization of 
the magnetic moments increases with decreasing V kd H 
and controls the extent to which rare earth impurities 
enhance damping. 

The orbit-orbit coupling (cf. Eq.|2J gives no contribu- 
tion for Gd +++ (4/ 7 ) as observed in the experiment [l2j |. 
As for the element Eu, it is believed from measurements 
of the paramagnetic susceptibilities that the ionic state is 
Eu++ (4/ 7 ) and not Eu+++ (4/ 6 ) [HIH. If this is the 
case then clearly this is a state with L/ = and it is the 
same as that of Gd +++ . Yb is also present in a double- 
ionized state [I]} and therefore doping with Yb ++ (4/ 14 ) 
should not increase damping. This result remains to be 
confirmed by experiment. For Eu there is an additional 



reason why its angular momentum is quenched. The first 
excited state of this latter element lies only about 400 K 
above the ground state and this can lift the degen- 
eracy of the ground state. The average orbital angular 
momentum will therefore be zero even though L 2 remains 
a good quantum number |28j . Hence our Hamiltonian 
from the outset reproduces the experimental results for 
Eu and Gd and predicts that doping with Yb should not 
change the damping. We next address the remaining 
rare earth elements. 

First, we outline the steps to derive the damping due 
to the orbit-orbit coupling term. We are only interested 
in the damping of the (i-moments of the transition metal, 
therefore it is advantageous to adopt a functional integral 
approach. Since our system is near equilibrium and far 
from the Curie point, we use the spin wave approxima- 
tion and expand the spin operators of the f-moments in 
terms of Boson operators where f ± ~Sj± iSJ. We 
keep only the first non-trivial terms. The integration of 
the conduction electrons is carried out exactly. After- 
ward we integrate the impurity variables, / and f\ also 
exactly but keep only quartic terms in d and d + . The re- 
maining effective action has now only the fields d and w 
and from their equations of motion the spin propagator 
(to - (r)m + (r')) of the d-moments, nv^ = ± iS%, can 
be determined. We use a Stratonovich-Hubbard transfor- 
mation to write this effective Lagrangian in terms of m . 
Then a stationary phase approximation of the functional 
generator allows us to determine the desired propagator 
and hence the damping. We finally compare the func- 
tional form of this result to that of LLG and discuss why 
the electronic (itinerant) picture of the host transition 
ions is essential. 

The fundamental quantity in our calculation is the gen- 
erating functional 

Z[ri*,r)] = Tre~f° dr f H -')*(- r ) m_ ( T )-'?( r ) m+ ( T )}. (5) 

where 77 and rj* are external sources and (3 is inverse 
temperature. The propagator, i.e. the connected two- 
point Green's function, of the volume mode of the 
transition metal ions is found by functional differenti- 
ations with respect to the external sources 77* and rj, 
(m+(T)m-{T')} c = 6 2 lnZ[rj*,ri}/6ri(T)Sri*{T'). It is cal- 
culated within a double random phase approximation 
(RPA2) method. The true single particle propagator of 
the d-bands is first found within a RPA in the presence 
of an effective field due to the conduction electrons and 
the impurities. In turn, the effect of the /-impurities on 
the conduction electrons is calculated within RPA. The 
resulting effective Lagrangian is now written in terms of 
m only 

C = --my£y W mjw - Tr In [G^ 1 + JCm] . (6) 

where G d ~\ai, a 2 ) = d T -e d + V 2 G c + tr k {G f G c BG c A} 
is the propagator of the d-electrons in the presence of 
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the conduction electrons and the rare earth impurity 
(o"j = 1,2 for spin up and spin down respectively). The 
quadratic term in m represents effective anisotropy and 
spin-charge interactions and is given by 



In the absence of impurities, these equations are to lowest 
order the time-dependent generalization of the Hartree- 



K. 



(7\ (T2CT3 £T4 

2t/4 



2Sl 



Fock equations derived by Anderson [2? 
RPA2 method, we solve for C(1221) 

( 7 )Ci22i(^;) = V"m-ii(w„)m 22 (^n + u>i) 



Using the 



—2 V Gf(G c BG c AG c ) rTlc72 Gf(G c BG c AG c ) a - 3rTi d rTll j i S l j 2Cr3 
-V i G c AG c G fG c BG c . 

Integrations over momentum and spin are implied in all 
these expressions. The different terms that appear in K 
are as follows: G c is the Green's function of the conduc- 
tion electrons in the mean field approximation 



G 



c 1 fk,cri,k , <t 2 ,t) = (d T +£ko-i - Mf) 
+i\(k,k')(2- gj)(J z f ) (k' x k y - k' y k x ^j <W 2) (8) 

which is off-diagonal in momentum due to the orbit-orbit 
coupling. Ek,o- now includes Zeeman terms due to the 
external field and the z-component of the field due to 
impurity. The propagator G/ is that of the /-ions in 
the presence of both the conduction electrons and the 
transition ions, GJ 1 (r) = d T + fj,fH + Tr kt „ {G C AG C B}. 
The A and B matrices are solely due to the presence of 
the impurity and represent the indirect coupling between 
the transition ions and the /-ions 

A (k', ay, k, (t 2 ) = B (k, <x x - k', a 2 )* = r <7+ CT2 - iA A+ (9) 



we have set Tq 



and A± , 
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k'xk 



k'k v 

-1), A = 
In the trace 



log term of the effective Lagrangian, the first nontrivial 
contribution is of order V 4 and is given by Fig. ^ The 
diagram with a single insertion of an /-propagator does 
not contribute due to the antisymmetry of the orbit-orbit 
coupling in the momentum space. Varying the effective 





FIG. 1: The first diagram that is contributing to the damping 
of the d-electrons due to the f-impurities through the conduc- 
tion electrons. 

action with respect to my gives four equations which 
can be averaged and differentiated with respect to the 
external sources to get the m propagators. We are only 
interested in C(1221) = (to 12 to 2 i) which is given by 



-(m 22 ) 



/Ciiy(m y )}C(1221) 
- /C 2 iyC(ij21)(m 22 ) - 



/Cn y C(zj21)(mi 2 )(10) 



n 

I 1 +H„, m ^2112(Wm)Wn(cJ„ +W m )TO 22 (w„ +L0 m + d>l) 

where u>i = (21 + 1)tt / (3 for integer I. If we ignore the im- 
purity interaction and replace the average values of the 
rriij by the Anderson solution, we recover the RPA result 
for the propagator of the magnetization. To include the 
impurities, we evaluate the d propagators, rriij, within 
RPA. In the low frequency limit, u << A << oj c , we 
find that the (retarded) propagator C R of the theory is 
proportional to (oj — luq + iau) . Here, A -1 is the life- 
time of the virtual d states , oj c denotes the frequency 
of the conduction electrons, and ujq is the ferromagnetic 
resonance frequency of the transition metal. This low 
frequency limit for the damping is similar to that of the 
LLG result 0|. The damping a in the spin-conserving 
channel is proportional to Jf(Jf + 1) ((gj — 2)\V\) 4 and 
is given by 



a = c|AF| 4 J / (J / + l)(2 
UAE 

x 



gj) 



(nraki 



(12 

Q(<»/) 



2 5 7r 3 (£- AE) 2 (E + AE) 2 18uj c 



Here n is the density of conduction electrons, c is the 
concentration of the /-impurities, and E±AE is the en- 
ergy of the up/down d states. These latter energies can 
be determined self-consistently as in the Anderson solu- 
tion |23 | and hence their form is not expected to depend 
strongly on the atomic number of the rare earth impurity 
at low concentrations. The explicit form of the function 
Q is not needed here but it represents contributions be- 
yond the 'mean' field approximation of the /-impurities 
and is given by Fig. ^ In Fig. |2 we show that the 
leading coefficient of the damping due to non-spin flip 
scattering (solid curve) is in very good agreement with 
the experimental results of Bailey et al. 12]. 

Finally we point out the reasons behind insisting on us- 
ing the itinerant electrons explicitly instead of the simpler 
s-d exchange interaction which accounts well for damping 
in permalloy [l6| . Using a localized-type Hamiltonian for 
the (i-moments 



(13) 



,)C(1221). 



instead of Eq. leads to a damping which differs signif- 
icantly from experiment (dashed curve in Fig. 0). This 
localized moment Hamiltonian however appears to de- 
scribe well damping in insulators such as heavy rare earth 
doped garnets 29]. In garnets, the hybridization cou- 
pling is smaller than in metals. Hence our result also 
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explains why the damping in rare-earth doped garnets is 
not as strong as in the rare-earth doped transition metals. 
The experimental measurements (triangles) clearly show 
that at room temperature non-spin flip scattering is more 
important than spin-flip scattering which only becomes 
important close to the critical temperature. Again, the 
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FIG. 2: Comparison of the normalized leading factor in the 
damping as a function of the rare earth impurity in Eq. Y\^\ 
(solid line) and Eq |13l (dashed line) to the data of Ref. Il2i 
The squares represent damping due to s-f coupling only, Eq. 
|5J without the orbit-orbit coupling. 

data is well reproduced by the orbit-orbit coupling and 
the relatively large increase in damping is due to the large 
virtual mixing parameter Vkd- In constrast, the s-f cou- 
pling (squares in Fig. |5J) is in conflict with experiment. 

In summary, we have shown that the damping in rare- 
earth doped transition metals is mainly due to an orbit- 
orbit coupling between the conduction electrons and the 
impurity ions. For near equilibrium conditions and in the 
low frequency regime this leads to damping for the uni- 
form mode that is of Gilbert form. The orbit-orbit mech- 
anism introduced here is much stronger than the spin- 
orbit based Elliott- Yafet-Kambersky mechanism since 
the charge-spin coupling at the host ion is of the order 
of 1-10 eV compared to 0.01 eV for spin-orbit coupling. 
The predicted increase of damping is proportional to V 4 
which in transition ions is of the same order as U the 
Coulomb potential. A localized model for the d-moments 
based on the s-d exchange is unable to account for the in- 
crease in damping in these doped systems as a function of 
the orbital moment of the rare-earth impurities. An ad- 
ditional test of this damping theory would be to measure 
the effect of a single rare earth element on the damping 
in various transition metals. Such experiments will pro- 
vide further insight into the dependence of damping on 
V and will improve our understanding of the itinerant 
versus localized pictures of magnetism. 
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